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ABSTRACT 

The  aerodynamic  motion  of  rectangular  plates  observed  in  18  free 
fall  drop  tests  is  analyzed  and  described.  In  most  of  the  drops,  the 
plates  autorotated;  the  autorotation  was  independent  of  the  plates' 
initial  attitude.  The  autorotation  generated  a  small  lifting  force 
which  caused  the  plates  to  experience  a  gliding  and  spiraling  motion. 
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SUMMARY 


The  aerodynamic  motion  of  rectangular  plates  observed  in  18  free- 
fall  drop  tests  is  analyzed  and  described.  In  most  of  the  drops,  the 
plate  autorotated;  the  autorotation  was  independent  of  the  plates'  ini¬ 
tial  attitude.  The  autorotation  generated  a  small  lifting  force  which 
caused  the  plates  to  experience  a  gliding  and  spiraling  motion. 

The  angular  rate  around  the  minimum  inertia  axis  of  each  plate 
was  obtained  from  camera  data.  A  curve  and  the  corresponding  equation 
were  used  in  fitting  the  data.  The  curves  obtained  indicate  that  the 
plates  exhibited  the  typical  characteristics  of  an  autorotating  body, 
with  a  peak  buildup  of  angular  rate  and  a  decay  to  a  steady-state  value 

Trajectory  data  were  obtained  for  all  but  two  of  the  tests.  Equa 
tions  were  derived  which  define  the  drag  coefficient  as  a  function  of 
velocity.  These  equations  indicate  a  decrease  in  drag  coefficient  with 
an  increase  in  velocity.  There  was  good  agreement  between  the  calcu¬ 
lated  velocities  and  the  experimentally  determined  velocities;  this  was 
an  additional  indication  that  the  plates  experienced  autorotational 
motion. 

The  maximum  dispersion  of  the  plates  which  were  dropped  from  2000 
feet  above  terrain  was  within  a  circle  1000  feet  in  radius;  the  plates 
dropped  from  10,000  feet  above  terrain  were  dispersed  within  a  circle 
6000  feet  in  radius.  In  general,  the  plates  suffered  little  damage  at 
impact,  and  most  of  them  impacted  on  their  longer  side. 
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FREE-FALL  ROTATION  AND  AERODYNAMIC  MOTION 
OF  RECTANGULAR  PLATES 


Introduction 

A  study  was  undertaken  to  investigate  the  free -flight  motion  of 
rectangular  plates  because  components  having  this  configuration  will 
reenter  the  earth's  atmosphere  as  parts  of  radioisotopic  power  genera¬ 
tors.  For  some  of  the  planned  missions,  it  is  desired  that  these  com¬ 
ponents  survive  the  reentry  environment  and  remain  intact  after  impact¬ 
ing  the  earth's  surface.  It  is  necessary  to  know  what  free-flight 
motion  the  plates  exhibit  during  the  terminal  phase  of  the  flight  in 
order  to  assess  the  plates  structural  integrity  and  trajectory  charac¬ 
teristics.  In  addition,  it  is  planned  to  investigate  the  free-flight 
motion  of  these  configurations  in  a  hypersonic  environment. 

The  first  part  of  the  study  was  a  drop-test  program  in  which 
rectangular  plates  were  freely  dropped  from  altitudes  of  2000  and 
10,000  feet  above  terrain.  These  tests  were  conducted  in  order  to  ob¬ 
tain  the  aerodynamic  motion  and  free-flight  trajectory  characteristics 
of  the  plates.  The  results  and  analysis  of  18  drop  tests  conducted  at 
the  Tonopah  Test  Range  in  Nevada  are  presented. 


General  Theory  and  Background 

This  study  was  undertaken  to  investigate  the  motion  exhibited  by 
a  rectangular  flat  plate  during  free  flight  and  to  determine  the  depend¬ 
ence  of  this  motion  on  initial  conditions.  Before  proceeding  into  the 
actual  study,  it  is  worthwhile  to  present  the  general  theory  which  de¬ 
scribes  the  free  motion  such  a  plate  should  undergo.  A  flat  rectangular 
plate  in  free  flight  will  generally  have  a  tumbling  motion.  Tumbling 
is  defined  in  this  report  as  a  free  body  rotation,  predominantly  about 
an  axis  perpendicular  to  the  airstream  but  not  always  because  a  combin¬ 
ation  of  motions  may  often  be  discerned.  However,  it  is  felt  that  a 
rectangular  flat  plate  would  generally  experience  an  autorotational 
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motion.  In  this  report,  autorotation  is  defined  as  rotation  about  a 
fixed  axis  which  is  perpendicular  to  the  airstream.  This  rotation 
increases  to  a  p^alc  value  and  then  achieves  a  steady  state  rate  without 
the  aid  of  an  external  force.  Reference  1  gives  an  excellent  treatment 
of  tumbling  motion  and,  therefore,  no  detailed  theories  or  descriptions 
of  this  motion  are  given  here.  It  is  believed  that  a  rectangular  plate 
will  continuously  autorotate  about  an  axis  which  is  free  in  flight. 

This  behaviour  may  be  compared  to  the  motion  of  a  strip  of  paper  of 
aspect  ratio  two  or  more .  When  dropped  in  still  air ,  the  strip  rotates 
and  descends  along  an  inclined  path.  Aerodynamically ,  this  motion  may 
be  explained  by  the  fact  that  turbulent  flow  requires  a  certain  time  to 
form.  This  is  illustrated  by  Figure  1. 


(b)  Actual 


Figure  1.  Comparison  of  potential  and  actual 
flow  over  a  strip  of  paper 
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Potential  flow  is  generally  formed  at  the  instant  that  flow  begins 
over  a  body  and  it  is  so  called  because  it  can  be  calculated  by  use  of 
the  potential  theory  of  aerodynamics.  In  Figure  1  (a),  the  steady-state 
potential  flow  around  the  strip  is  shown  in  five  consecutive  positions. 
In  position  (a)l  no  torque  is  exerted  by  the  flow  on  the  strip.  In 
position  (a) 2,  the  theoretical  flow  has  two  stagnation  points,  S,  of  no 
velocity  and  of  maximum  pressure.  The  position  of  the  strip  is  such 
that  the  stream  exerts  a  clockwise  driving  torque  on  the  strip.  The 
flow  in  position  (a)3  is  symmetrical,  and  no  torque  is  exerted.  Posi¬ 
tion  (a)4  appears  the  same  as  position  2(a),  but  the  torque  is  now  act¬ 
ing  counterclockwise  and,  since  the  strip  is  still  rotating  clockwise, 
amounts  to  a  retarding  torque.  In  the  quarter  revolution  between  (a)l 
and  (a) 3,  the  flow  performs  positive  work  on  the  strip,  but  between  (a) 3 
and  (a) 5,  it  performs  exactly  the  same  amount  of  negative  work  on  it. 

For  potential  flow  where  positions  (a)2  and  (a)4  are  identical, 
the  stream  performs  zero  work  during  a  full  revolution,  and  neither 
damping  nor  self-excitation  takes  place.  However,  the  actual  flow  dif¬ 
fers  from  the  potential  flow  and  is  shown  in  Figure  1(b).  Between 
(b)l  and  (b)3,  turbulence  behind  the  plate  builds  up  gradually,  and 
between  (b)3  and  (b)5  this  turbulent  air  is  washed  away  with  the  stream. 
Because  of  this,  positions  (b)2  and  (b)4  differ;  (b)2  resembles  the 
potential-flow  position  much  more  than  does  (b)4.  The  effect  of  tur¬ 
bulence  on  the  plates  in  position  2  is  primarily  to  diminish  the  torque 
exerted  by  the  stream  on  the  plate.  Similarly,  in  position  (b)4  the 
retarding  torque  is  less  than  that  of  (a)4,  but  the  "less"  is  more  pro¬ 
nounced  in  position  4  than  in  position  2.  Therefore,  there  is  more 
positive  work  performed  between  (b)l  and  (b)3  than  there  is  negative 
work  between  (b)3  and  (b)5,  and  net  positive  work  results  from  a  full 
revolution.  This  gives  rise  to  the  self -excited  rotation  or  "autorota¬ 
tion"  of  the  plate. 

Besides  a  torque,  the  flow  also  exerts  a  sidewise  force  on  the 
plate.  The  flow  in  Figure  1  (b)2  resembles  the  flow  around  an  airplane 
wing,  and  the  plate  experiences  a  lift  force  to  the  right.  Similarly, 
in  (b)4  there  is  a  lift  force  to  the  left,  but  this  force  is  smaller 
than  (b)2  because  the  turbulence  is  much  further  developed  ("stalled" 
airplane  wing).  Thus,  there  is  a  net  lift  force  to  the  right  during 
a  full  revolution.  This  explains  the  inclined  fall  of  the  paper  strip. 
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The  preceding  discussion  has  provided  a  qualitative  view  of  the 
existence  of  autorotation.  This  discussion  is  not  complete  because 
other  forces  can  contribute  to  autorotation.  This  introduction  into 
the  autorotation  phenomenon  was  presented  so  that  the  data  presented 
in  this  report  will  be  more  easily  understandable  in  terms  of  autoro¬ 
tation.  The  plates  in  this  study  have  six-degrees-of-freedom  available, 
and,  consequently,  their  motion  will  not  be  purely  autorotational  but 
will  be  influenced  by  several  possible  impressed  aerodynamic  forces 
which  can  induce  a  complex  motion,  particularly  when  gyroscopic  effects 
are  added. 


Description  of  Tests 

A  total  of  18  free-fall  drop  tests  was  conducted  in  which  rectan¬ 
gular  flat  plates  were  freely  dropped  in  various  attitudes.  These  tests 
were  conducted  from  altitudes  of  2000  and  10,000  feet  above  terrain. 

A  Kaman  H-43  helicopter  was  used  as  the  drop  helicopter,  and  a  Bell  G3B 
helicopter  was  used  for  photographic  coverage.  During  the  test,  the 
H-43  helicopter  was  flown  at  approximately  zero  air  relative  velocity, 
and  the  plates  were  dropped  from  the  rear  of  the  helicopter .  This  was 
done  in  order  to  minimize  any  initial  disturbance  or  motions  that 
might  have  been  imparted  to  the  plates  by  the  downwash  or  wake  of  the 
helicopter.  As  a  result  of  the  precautions  described  above,  it  is 
believed  that  the  initial  motion  of  the  plates  were  minimally  influenced 
by  the  helicopter's  wake  or  motion. 

The  tests  were  conducted  at  Sandia  Laboratories'  Tonopah  Test 
Range,  Tonopah,  Nevada,  which  is  quite  heavily  instrumented  with  ground 
cameras  and  radars.  Both  radar  and  ground  optical  tracking  cameras  were 
used  to  obtain  trajectory  and  photographic  data.  In  addition,  the  chase 
helicopter  had  a  photographer  on  board  who  hand-tracked  the  drops  with 
a  16 -mm  movie  camera. 

Two  different  rectangular  plates  were  used:  a  32-  x  27-  x  0.69- 
inch  plate  with  hemicylindrical  edges  and  a  48-  x  18-  x  0.50-inch  plate 
with  square  edges.  The  physical  characteristics,  weights,  and  moments 
of  inertia  are  shown  in  Table  I . 
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TABLE  I 

Physical  Characteristics  of  Plates 


Moments  of  Inertia 

Dimensions 

Weight 

Ixx 

Iyy 

Izz 

(inches) 

(lbs) 

(slug-ft2) 

(slug-ft2) 

(slug-ft2) 

No.  1  -  32  x  27  x  0.69 

40±l/4 

1.199 

0.5098 

0.693 

No.  2  -  48  x  18  x  0.50 

44±l/4 

2.08 

0.257 

1.812 

The  test  program  was  conducted  in  two  separate  series.  In  the 
first  series,  eight  drop  tests  were  conducted  with  Plate  No.  1  on  July 
11  and  12,  1966.  The  plate  was  tricolored  (black,  yellow,  and  orange) 
to  make  it  easily  visible  and  to  aid  in  obtaining  angular  rates  (Figure 
2).  The  second  series  was  conducted  on  August  15  and  16,  1966;  six 
drops  were  conducted  with  Plate  No.  1  and  four  with  Plate  No.  2.  Both 
plates  had  varicolored  geometric  designs  to  aid  in  tracking  and  obtain¬ 
ing  angular  rates  (Figure  3).  The  plates  were  dropped  2000  and  10,000 
feet  and  at  seven  different  initial  attitudes  (Table  II). 


Test  Analysis  and  Discussion 


Aerodynamic  Motion 

The  main  objectives  of  the  test  were  to  investigate  the  type  of 
motion  that  a  rectangular  plate  experiences  in  free  flight  and  to  de¬ 
termine  if  this  motion  is  independent  of  the  initial  conditions  or 
disturbances.  Of  secondary  importance  were  the  determination  of  the 
impact  characteristics  of  a  plate  and  its  terminal  velocity.  The  tests 
were  conducted  at  Sandia  Laboratories  Tonopah  Test  Range  in  Nevada. 
Instrumentation  consisted  of  optical  range  cameras  for  trajectory  and 
photographic  coverage;  in  addition,  photographic  data  were  provided  by 
a  chase  helicopter.  The  motion  on  all  18  drop  tests  was  analyzed  by 
studying  the  film  data.  The  plates  were  released  in  seven  different 
attitudes,  with  some  of  the  attitudes  repeated.  The  two  altitudes  of 
release  were  2000  and  10,000  feet  above  terrain. 

In  each  test  except  Test  No.  258-19,  the  plates  went  into  an 
autorotational  mode  (see  Table  II).  In  Tests  258-13,  -16,  and  -17, 
the  plates  were  dropped  in  a  flat  attitude  and  remained  in  this  attitude 
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TABLE  II 


Test  Summary 


Test  No. 
258-10 


Initial  Attitude 


I 


Altitude 
Above  Terrain 


mi 

2000 


Final  Motion 
Aut or ota ting 


258-11  Edge  on 


2000  Autorotating 


258-12 


Autorotating 


258-13  Flat 


10,000 


1 


Flat  till  90 
seconds,  then 
started  auto¬ 
rotating 


(continued  on  next  page) 
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TABLE  II  (cont) 


Test  No. 
258-14 


258-15 


258-16 


_ Initial  Attitude 

Edge  on 


Altitude 
Above  Terrain 
(ft) _ 

10,000 


Rotate  plate  about  axis  xx 


10,000 


Flat  10,000 


v 


Final  Motion 
Autorotating 


Autorotating 


Flat  till  50 
seconds,  then 
started  auto¬ 
rotating 


258-17  Flat 


Flat  till  55 
seconds,  then 
started  auto¬ 
rotating 


(continued  on  next  page) 
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TABLE  II  (cont) 


Test  No. 
258-18 


258-19 


258-20 


258-21 


_ Initial  Attitude 

Flat 


32" 


Altitude 
Above  Terrain 


(ft) 

2000 


Final  Motion 
Autorotating 


Edge  on 


Flat 

(broadside) 
till  impact 


Rotate  plate  about  axis  xx 


2000  Autorotating 


(continued  on  next  page) 
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TABLE  II  (cont) 


Test  No. 
258-22 


258-23 


258-24 


258-25 


_ Initial  Attitude 

Edge  on 


Altitude 
Above  Terrain 


2000 


Final  Motion 
Autorotating 


I 


Rotation  about  yy  axis 


2000  Autorotating 


Edge  on 


Horizontal  sliding-out  motion 
(long  edge) 


2000 


10,000 


(continued  on  next  page) 


Autorotating 


Autorotating 
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TABLE  II  (cont) 


Test  No. 
258-26 


258-27 


Initial  Attitude 


Horizontal  sliding-out  motion 
(short  edge) 


Altitude 
Above  Terrain 

(ft) 

10,000 


Rotation  about  yy  axis 


10,000 


Final  Motion 
Autorotating 


Autorotating 


for  approximately  50  to  90  seconds  before  they  went  into  an  autorota- 
tional  mode.  The  flat  attitude  (or  a  =  90°)  is  a  statically  stable 

■jt 

position  for  the  plates.  In  each  case  in  which  a  plate  autorotated, 
the  plates  entered  the  autorotating  mode  by  building  up  a  rocking  am¬ 
plitude  in  a  vertical  plane.  The  plate's  first  full  revolution  of 
rotation  was  completed  by  rotating  in  a  plane  inclined  45°  to  60°  from 
the  vertical.  Simply,  then,  the  plate  did  not  have  enough  kinetic 
energy  to  rotate  over  "top  dead  center."  However,  rotation  was  achieved 
about  an  axis  inclined  to  a  position  requiring  less  potential  energy  to 
support  the  rotation.  After  rotation  was  well  established,  the  plate 
in  some  cases  returned  to  vertical. 


There  were  cases  where  the  plates  started  autorotating  almost 
immediately.  This  occurred  mostly  when  an  initial  angular  rate  was 
imparted  to  the  plate . 
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The  motion  of  the  plates  was  irregular.  This  is  understandable 
since  Plate  No.  1  was  nearly  square  and  had  a  low  aspect  ratio.  A  more 
stable  rotation  can  be  achieved  by  increasing  the  aspect  ratio,  i.e., 
making  a  longer  rectangle.  This  was  the  case  for  Plate  No.  2,  whose 
aspect  ratio  was  2.67  compared  to  Plate  No.  l's  aspect  ratio  of  1.12. 
For  Plate  No.  2,  the  rotational  rate  was  higher  and  more  nearly  regular 
than  it  was  for  Plate  No.  1.  This  regular  and  irregular  motion  is  a 
function  of  the  dynamic  moment  which,  in  turn,  depends  on  the  moment  of 
inertia.  This  is  illustrated  by  the  dynamic -moment  equation 


where : 


1a  ~  %)  +  M(A) 
a  =  Angle  of  attack 

=  Static  aerodynamic  moment,  for  ot  constant 

=  Net  driving  or  damping  moment  resulting  from 
'A7  a  plate's  rotation;  will  be  referred  to  as  auto- 
rotational  moment . 

I  =  Moment  of  inertia  about  axis  of  rotation. 


The  above  equation  can  be  analyzed  qualitatively  considering  Figure  4. 


Figure  4.  Variation  of  aerodynamic  and 
autorotational  moments  with 
angle  of  attack 
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The  above  curve  shows  the  variation  of  the  moments  and  with 

a.  If  the  integral  f  M(fl)dahas  a  positive  value,  it  represents  the 
motive  energy  of  the  moment  during  the  half  revolution,  which  means 

that  the  net  autorotational  moment  is  driving,  i.e.,  =  M(A)p  -  M(a)2 

where  now  M(A)p  >  M(A)2*'  For  steady  state  motion 


(2) 


that  is  the  motive  energy  M^a)-^  is  now  balanced  by  the  resisting  energy 
M(a)2*  These  moments  are  represented  by  the  shaded  areas  in  Figure  4. 

In  the  motive  interval,  the  autorotational  moment  is  causing  autorota¬ 
tion,  and  for  the  rest  of  the  period  it  is  retarding  autorotation.  In 
the  motive  interval,  a  plate  must  store  up  enough  kinetic  energy  to 
clear  the  second  interval.  Consequently,  the  motion  cannot  be  sustained 
unless  the  moment  of  inertia  of  the  plate  is  sufficiently  high;  other¬ 
wise,  the  rotation  stops  in  a  position  where  the  plate  starts  to  slide. 
In  general,  this  explains  why  the  plates  with  an  aspect  ratio  near  1 
have  an  irregular  steady  motion,  while  the  plates  with  an  aspect  ratio 
of  2.67  have  a  more  regular  rotation;  this  results  from  the  reduced 
chord.  In  general,  it  can  be  said  that  the  basic  free -flight  motion  of 
the  plates  was  an  autorotational  mode,  with  the  possibility  of  a  broad¬ 
side  mode.  The  autorotational  mode  is  much  more  predominant;  the  broad¬ 
side  mode  is  a  special  case.  Analysis  of  the  photographic  data  clearly 
indicates  that  the  motion  of  the  plates  is  independent  of  the  initial 
conditions  and  that  the  final  motion  of  a  plate  is  either  autorotational 
or  a  broadside  attitude  (see  Table  II) . 


The  autorotational  rate  of  the  plates  was  read  from  the  photo¬ 
graphic  film.  It  was  difficult  to  obtain  a  single-degree-of -freedom 
rotational  rate  from  the  film  because  of  the  compound  motion.  This  com¬ 
pound  motion  consists  of  the  plate  rotating  about  a  tilt  axis  and,  at 
the  same  time,  precessing  and  translating  (Figure  5).  The  plate  is 
rotating  about  the  yy  axis  which  is  perpendicular  to  the  zz  axis;  the 
zz  axis  is  tilted  from  a  vertical  axis.  Also,  there  is  a  precessional 


or 


/C  • 

NOTE:  This  autorotational  net  moment  can  be  either  damping  or 
driving,  i.e.,  £  M(A)2‘ 
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motion  superimposed;  this  motion  has  an  angular  rate,  Q,  around  the 
vertical  axis  . 


Figure  5 .  Diagram  of  plate’s  compound  motion 

Average  angular  rotational  rates  around  the  plate's  axis  of  mini¬ 
mum  inertia  were  obtained  from  the  film  data  as  a  function  of  time  for 
the  plates;  the  rates  are  shown  in  Figure  6  for  some  of  the  tests. 

Since  these  rates  are  fairly  well  scattered,  a  least  square  computer 
technique  was  used  to  fit  the  angular  rate  data.  Two  methods  of  solu¬ 
tion  were  used,  "normal"  least  squares  and  orthogonal  polynomials;  the 
normal  least  square  method  gave  the  best  fit.  (See  Appendix  for  tech¬ 
nique)  .  Using  this  technique,  angular  rate  curves  versus  time  were 
obtained  for  two  of  the  tests  and  are  shown  in  Figures  7  and  8  along 
with  their  representative  equation.  These  curves  exhibit  the  typical 
characteristics  of  an  autorotating  body  with  a  peak  buildup  of  angular 
rate  and  a  decay  to  a  steady-state  value  of  approximately  5  rps .  In 
general,  the  dynamic  data  obtained  from  these  tests  are  fairly  qualita¬ 
tive,  so  no  specific  results  or  analysis  can  be  obtained.  However,  they 
do  show  that  the  plates  have  a  compound  autorotational  motion  and  that 
the  final  motion  of  a  plate  is  independent  of  its  initial  condition. 
These  conclusions  are  based  on  analysis  of  the  film  data.  In  order  to 
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Plate  s  average  angular  rate 
obtained  from  camera  data 


Figure  7 .  Curve  fit  of  angular  rate  data 
for  one  drop  test 


0  5  10  15  20  25  30 

Time  (sec) 

Figure  8.  Curve  fit  of  angular  rate  data 
for  one  drop  test 

do  a  complete,  detailed  analysis  of  this  motion,  the  full  6-degree-of - 
freedom  aerodynamic  coefficients  are  required.  This  is  planned  for 
future  studies  where  the  dynamic  aerodynamic  coefficients  will  be  ob¬ 
tained  from  wind-tunnel  tests  and  then  used  to  predict  analytically  the 
motion  of  the  plate. 


Drag  Coefficients 

The  drag  coefficient  (Cp)  of  the  plates  was  calculated  as  a  func¬ 
tion  of  time,  velocity,  and  Mach  number  from  the  trajectory  data*  ob¬ 
tained  from  the  ground-based  tracking  cameras .  The  referenced  area 
used  for  the  coefficient  was  based  on  the  6  ft^  broadside  area  of  the 
plate.  A  typical  drag  coefficient  versus  time  plot  for  three  drop 
tests  is  shown  in  Figure  9.  The  plate  in  Test  258-13  remained  in  a 
flat,  rocking  attitude  for  90  seconds  before  autorotating ;  in  Tests 
258-16  and  258-17,  the  plates  remained  in  a  flat,  rocking  attitude  for 
roughly  50  seconds  before  autototating .  As  indicated  in  Figure  9,  the 
drag  remains  fairly  constant  with  time  until  the  plates  go  into  an 
autorotating  mode,  at  which  time  the  curve  shows  oscillations  in  the 
data.  These  oscillations  result  from  the  variations  in  the  trajectory 
caused  by  the  gliding  and  spiraling  effects  of  the  plate.  Also,  some 
of  these  variations  result  from  the  manner  in  \4iich  the  data  are 
smoothed  and  reduced  in  order  to  obtain  the  drag  coefficient. 

The  drag-coefficient  data  were  used  to  obtain  an  equation  for 
Cd  as  a  function  of  velocity.  The  drag  data  for  each  test  were  curve - 
fitted,  using  a  least-squares  fit  computer  program  technique.  This 
method  is  presented  in  the  Appendix.  First  and  second  order  equations 
were  obtained  from  this  curve-fit  program.  In  general,  the  first-order 
equation  was  satisfactory,  and  the  equations  were  of  the  form 

CD  =  Kx  -  K2V  ,  (3) 

where 

V  =  velocity 

Kl>  ^2  are  arbitrarY  constants. 

A  typical  drag  coefficient  versus  velocity  curve  and  its  fit  are  shown 
in  Figure  10.  For  each  test,  and  K/?  had  different  values;  however, 
for  the  tests  considered,  the  constants  ranged  in  value  as  shown  below: 

1.50  <  KX  <  4.00 

0.0050  <  K2  <  0.0300. 


The  complete  trajectory  data  (i.e.,  range,  altitude,  velocity, 
acceleration,  etc.,  versus  time)  were  obtained  for  all  tests  except 
Tests  258-19  and  258-22.  No  data  were  obtained  for  these  tests  because 
the  plates  were  mis tracked. 
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Figure  9.  Variation  of  drag  coefficient 
of  plate  during  its  flight 


An  attempt  was  made  to  determine  if  the  constants  in  Equation  (3) 
were  dependent  in  any  way  on  the  initial  attitude  of  the  plates  and  if 
the  constants  could  be  correlated  from  test  to  test.  No  dependency  or 
correlation  was  found.  The  drag  coefficient  was  calculated  for  a 
velocity  of  90  fps  for  each  initial  attitude.  These  average  drag  val¬ 
ues  are  shown  below  as  a  function  of  initial  attitude. 


-Drag  coefficient 


Initial 

Attitude 

CD 

Flat 

0.980 

Rotating 

0.940 

Sliding  out 

1.025 

Edge  on 

0.955 

These  drag  coefficients  compare  quite  well  with  the  drag  coefficient  of 
a  plate  normal  to  the  flow  at  subsonic  speeds  (see  Reference  2). 

Equation  (3)  and  Figure  10  show  the  drag  decreasing  with  increas¬ 
ing  velocity.  At  first,  this  might  appear  aerodynamically  perverse 
since,  generally,  drag  increases  with  increasing  velocity.  However,  for 
a  rotating  flat  plate  for  a  given  velocity  or  Reynolds  number,  the  drag 
will  decrease  with  increasing  velocity.  It  is  important  to  point  out 
that  this  is  true  only  for  a  certain  limited  velocity  range  and  that 
Equation  (3)  applies  only  for  a  certain  velocity  range.  In  this  case 
the  decrease  of  drag  coefficient  with  velocity  can  be  explained  by 
analyzing  the  local  flow  over  the  rotating  plate.  This  flow  is  analo¬ 
gous,  in  the  limit  (i.e.,  a  large  rotational  rate,  u ,  to  the  flow  over 
a  cylinder.  Therefore,  for  certain  w/V  velocity  ratios,  the  separation 
on  the  plate  is  delayed,  and,  consequentially,  the  wake  (pressure)  drag 
decreases . 


Figure  10.  Variation  of  plate's  drag  coefficient 
with  velocity 
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Equations  of  the  form  of  Equation  (3)  were  obtained  and  plotted 
for  each  drop  test  and  are  shown  in  the  appendix. 

Velocity 

Velocity  data  were  obtained  for  each  drop  test  except  Tests  258-19 
and  258-22.  These  data  are  derived  from  the  trajectory  data  provided 
by  the  tracking  cameras.  Figures  11  and  12  show  some  typical  variations 


Time  (sec) 

Figure  11.  Typical  variation  of  plate  velocity  with  time 
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Figure  12.  Typical  variation  of  plate  velocity  with  time 


of  velocity  with  time.  In  Figure  12,  the  oscillations  in  the  velocity 
versus  time  curve  result  from  the  smoothing  and  data  techniques  used 
in  obtaining  the  velocity. 

Prior  to  the  actual  tests,  the  terminal  velocity  for  the  plate 
was  calculated  as  a  function  of  altitude,  this  was  done  for  various 
plate  attitudes  (Figures  13  and  14) .  In  addition,  Figure  13  gives  the 
velocity  for  the  10,000-foot  drops,  and  Figure  14,  gives  the  veloci¬ 
ties  for  the  2000-foot  drops.  It  is  evident  from  Figures  13  and  14 
that  the  velocities  reached  in  the  tests  fall  between  the  two  curves 
which  indicate  plate  rotation  about  its  short  and  long  axes  (W/CjjA  = 
6.06  and  9.5);  this  is  where  the  data  should  fall  since  the  actual 
rotational  motion  of  the  plates  was  composed  of  these  two  motions . 
However,  there  are  two  exceptions  in  Figure  14  (2000-foot  drops);  the 
velocities  from  Tests  18  and  23  do  not  fall  between  these  two  curves; 
this  could  be  the  result  of  an  error  in  the  trajectory  data. 

In  general,  the  terminal  velocities  obtained  from  the  tests  agree 
quite  well  with  the  predicted  velocities  as  shown  in  Figures  13  and  14. 

Altitude /Range  Data 

One  of  the  items  of  interest  in  the  drop  tests  of  the  plates  was 
the  range /altitude  data.  These  data  are  helpful  in  determining  the 
dispersion  of  the  plates.  Since  the  plates  are  autorotating ,  circula¬ 
tion  is  induced  which  generates  a  lift  force.  This  force  tends  to 
impart  a  gliding  or  spiraling  trajectory  to  the  plates.  If  a  gliding 
trajectory  is  imparted  to  a  plate,  it  will  generally  sail  and  have  a 
greater  range  than  if  it  were  falling  without  a  lifting  influence.  The 
altitude  versus  range  data  were  plotted  for  all  the  tests,  (Figures 
15  and  16),  and  some  of  the  spiraling  trajectories  are  evident  in  these 
figures.  It  also  can  be  seen  that  the  maximum  range  for  the  10,000- 
foot  drops  was  6000  feet,  and  for  the  2000-foot  drops,  it  was  1000  feet. 
Using  the  altitude /range  data,  it  was  determined  that  for  the  10,000- 
foot  drops  the  maximum  L/D  developed  by  the  plate  was  0.70  and  the 
minimum  was  0.333.  For  the  2000-foot  drops,  the  maximum  L/D  of  the 
plate  was  0.400  and  the  minimum  was  0.100.  Essentially,  these  L/D 
ratios  show  that  some  appreciable  lift  was  generated  by  the  rotating 


“X - 

"Except  Nos.  19  and  22;  no  data  were  obtained  for  these  drops. 
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Figure  14.  Terminal  velocity  for  plates  dropped  from 
2000  feet  above  terrain 
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flat  plates;  since  the  autorotational  motion  was  irregular,  these  L/D 
ratios  were  fairly  low. 

An  attempt  was  made  to  see  if  there  was  a  correlation  of  range 
(maximum  or  minimum)  with  initial  attitude;  as  was  expected,  none  was 
obtained.  The  impact  range  for  the  various  initial  attitudes  is  shown 
in  Figure  17 . 


Impact  range  (ft  x  10"^) 


Figure  17 .  Impact  range  for  various  initial  attitudes  of  plates 
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Critical  Reynolds  Number 

It  is  worthwhile  to  discuss  the  critical  Reynolds  number,  and 
to  determine  if  the  plates  in  these  tests  were  in  turbulent  or  laminar 
flow  (below  or  above  critical  Reynolds  number) .  It  is  generally  as¬ 
sumed  that,  for  steady-state  flow  on  a  flat  plate  at  a  =  0°  with  care¬ 
fully  controlled  conditions  and  a  very  smooth  surface,  the  critical 
Reynolds  number  is  generally  5  x  10^  .  For  the  same  plate  at  a  =  90°, 
the  flow  is  generally  considered  to  be  turbulent  and  the  critical 
Reynolds  number  is  not  considered.  The  cases  of  a  =  0°  and  90°  are 
the  two  extremes  for  a  tumbling  body,  and  the  flow  on  a  tumbling  body 
is  not  in  steady  state.  In  addition  to  the  stream  turbulence,  the 
transition  Reynolds  number  depends  on  the  manner  in  which  the  pressure 
in  the  flow  and  the  velocity  vary  along  the  surface.  When  the  pressure 
gradient  is  favorable  (positive),  the  critical  Reynolds  number  increases. 
A  flow  deceleration  along  the  surface  decreases  the  Reynolds  number  for 
transition.  In  the  case  of  the  tumbling  or  autorotating  flat  plate, 
the  flow  on  either  side  of  the  plate  reverses  direction  twice  per  revo¬ 
lution;  hence,  the  critical  Reynolds  number  should  be  lower  than  5  x  10^ 
for  the  dynamic  case.  It  is  assumed,  therefore,  that  the  flow  on  the 
tumbling  flat  plate  is  predominantly  turbulent,  and  therefore  the  drag 
coefficient  is  for  turbulent  flow.  Even  under  steady-flight  conditions 
and  at  ambient  temperature  the  plate  would  be  in  turbulent  flow  (See 
Figure  18).  Based  on  the  plate 1  s  dimension  of  27  inches,  the  Reynolds 
number  would  be  greater  than  500,000. 

Impact  Characteristics 

An  item  of  interest  and  concern  in  these  tests  was  the  condition 
of  Plate  No.  1  at  impact.  Other  items  of  interest  were  the  attitude 
and  velocity  of  the  plates  at  impact.  The  impact  velocities  were  ob¬ 
tained  from  the  trajectory  data  and  are  shown  in  Figure  19.  The  veloc¬ 
ities  range  from  a  minimum  of  73  ft /sec  to  a  maximum  of  88  ft/sec. 

The  average  impact  velocity  was  83  fps .  As  can  be  seen  in  Figure  19, 
the  impact  velocity  was  fairly  independent  of  the  release  altitude. 

The  kinetic  energy  of  the  plate  at  impact,  based  on  an  impact  velocity 
of  88  fps,  was  4810  ft-lbs.  The  TNT  equivalent  of  this  kinetic  energy 
is  0.05  ounces.  Obviously,  the  plate's  energy  at  impact  is  quite 
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Figure  18.  Reynolds  number  variation 
with  altitude  for  two  drop 
tests 
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Figure  19.  Plate  impact  velocity  for  the  various  drops 
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small,  and  damage  to  a  structurally  designed  plate  should  be  either 
small  or  nonexistent. 

Figure  20  shows  two  plates  as  they  appeared  after  impact.  Table 
III  gives  the  damage  to  and  attitude  of  each  plate  at  impact.  Most  of 
the  plates  impacted  on  the  long  side;  this  indicates  that  the  free- 
flight  motion  of  the  plate  was  predominantly  rotation  about  the  mini¬ 
mum  axis  of  rotation  which  in  this  case  is  the  xx  axis  (Table  III) . 

The  indentations  made  by  the  plates  were  generally  from  3-  to  5-inches 
deep  and  approximately  4-inches  wide.  The  plates  in  general  did  not 
sustain  severe  or  extensive  damage;  in  fact,  for  Tests  258-16  and  -17, 
the  plates  used  were  the  same  ones  originally  tested  in  Tests  258-13 
and  -14.  In  general,  it  is  felt  that  the  impact  loading  on  the  plates 
was  not  severe,  and  that  a  properly  designed  structural  plate  would 
encounter  little,  if  any,  damage. 


Summary  and  Conclusions 

The  results  of  these  tests  show  that  the  free-flight  motion  of 
these  rectangular  plates  displayed  compounded  autorotational  charac¬ 
teristics.  This  motion  was  the  same  (in  general)  for  each  drop,  and 
was  independent  of  the  initial  conditions .  The  predominant  rotation 
of  the  plate  was  about  its  minimum  axis  of  inertia;  this  motion  was 
further  compounded  by  the  rotation  about  another  vertical  axis,  thus 
introducing  gyroscopic  effects.  The  motion  data  was  obtained  primarily 
from  photographic  data.  Angular  rates  were  obtained  from  the  photo¬ 
graphic  data  and  fitted  with  a  least  square  curve  which  gave  the  angu¬ 
lar  rate  as  a  function  of  time.  This  curve  exhibited  the  typical  auto¬ 
rotational  rate  characteristics,  with  a  peak  build-up  in  angular  rate 
followed  by  decay  to  a  steady-state  value.  For  the  curves  obtained, 
this  steady-state  rate  was  5  rps . 

Trajectory  data  were  obtained  for  all  tests  except  two.  By  using 
these  data,  drag  coefficients  for  the  plates  were  obtained.  Also  deter¬ 
mined  were  expressions  which  related  the  drag  coefficient  as  a  function 
of  velocity.  These  expressions  showed  a  decrease  of  drag  coefficient 
with  velocity.  This  was  attributed  to  the  reduction  of  turbulence, 
i.e.,  the  reduction  of  separation  on  the  plate,  which  in  turn  reduced 
the  pressure  drag.  These  expressions  were  obtained  for  each  drop  test. 


Figure  20 .  Plates  after  impact 


Figure  20.  Plates  after  impact 


TABLE  III 


Impact  Characteristics 


Condition  of  Plate 


Test  No. 

Configuration 

Attitude 

at  Impact 

Following  Impact 

258-10 

Plate 

No.  1 

Landed 

on 

No.  1  side 

Slightly  damaged 

258-11 

Plate 

No.  1 

Landed 

on 

corner 

Plate  split 

258-12 

Plate 

No.  1 

Landed 

on 

No.  2  side 

Minor  damage 

258-13 

Plate 

No.  1 

Landed 

on 

No.  1  side 

Good 

258-14 

Plate 

No.  1 

Landed 

on 

No.  2  side 

Good 

258-15 

Plate 

No.  1 

Landed 

on 

No.  2  side 

Fair 

258-16 

Plate 

No.  1 

Landed 

on 

No.  2  side 

Fair 

258-17 

Plate 

No.  1 

Landed 

on 

No.  2  side 

Destroyed 

258-18 

Plate 

No.  1 

Landed 

on 

No.  2  side 

Good 

258-19 

Plate 

No.  1 

Landed 

on 

corner 

Badly  bowed 

258-20 

Plate 

No.  1 

Landed 

on 

corner 

Badly  bowed  and 
splintered 

258-21 

Plate 

No.  2 

Landed 

on 

No.  2  side 

Not  applicable 

258-22 

Plate 

No.  2 

Landed  on  corner  and 
No.  1  side 

Not  applicable 

258-23 

Plate 

No.  2 

Landed 

on 

No.  2  side 

Not  applicable 

258-24 

Plate 

No.  2 

Landed 

on 

No.  1  side 

Not  applicable 

258-25 

Plate 

No.  1 

Landed  on  corner  and 
No.  1  side 

Badly  damaged 

258-26 

Plate 

No.  1 

Landed 

on 

No.  2  side 

Badly  damaged 

258-27 

Plate 

No.  1 

Landed 

on 

No.  2  side 

Good 

No.  1  side 


No.  1  side 


Plate  No.  1 


Plate  No.  2 
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The  velocity  for  each  drop  test  was  also  obtained.  The  velocity 
as  a  function  of  altitude  obtained  from  the  tests  was  compared  against 
that  predicted,  and  the  agreement  was  good.  This  comparison  further 
showed  that  the  plates  exhibited  autorotational  motion. 

Range  and  altitude  data  obtained  from  these  tests  provided  infor¬ 
mation  on  dispersion  characteristics  of  the  plates  .  The  maximum  dis¬ 
persion  of  the  plates  which  were  dropped  from  2000  feet  above  terrain 
was  within  a  circle  1000  feet  in  radius;  the  plates  dropped  from  10,000 
feet  above  terrain  were  dispersed  within  a  circle  6000  feet  in  radius. 
These  data  also  showed  that  the  autorotational  motion  of  the  plate 
generated  some  appreciable  lift  forces. 

Based  on  the  Reynold's  number  criterion,  it  was  determined  that 
the  plates  were  in  turbulent  flow  throughout  the  flight. 

The  impact  data  on  the  plates  showed  that  the  model  plates  gen¬ 
erally  suffered  little  at  impact.  Impact  velocities  for  the  plates 
were  obtained,  and  it  was  determined  that  the  plates  impacted  predom¬ 
inantly  on  the  long  side. 
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Appendix 


Methods  Used  in  Obtaining  Equations  of  Angular  Rate 
Versus  Time  and  of  Drag  Coefficient  as  a 
Function  of  Velocity 

This  appendix  presents  the  method  used  to  obtain  equations  of 
angular  rate  versus  time  and  of  drag  coefficient  as  a  function  of  veloc¬ 
ity.  The  first  order  equations  of  drag  coefficient  versus  velocity 
obtained  for  each  test  are  presented  in  Table  A-I,  and  plots  are  shown 
in  Figures  A-l  through  A-32.  Shown  in  these  plots  are  the  linear-fit 
line  of  the  data  and  the  error  that  exists  between  the  computed  value 
and  the  actual  data  value.  As  is  evident,  the  error  is  generally  small. 

The  following  is  a  description  of  the  method  used.  This  method 
was  formulated  by  the  Mathematical  Computing .Services  Division  of  Sandia 
Laboratories ,  Albuquerque . 

Curve -Fit  Technique 

A  least-squares  technique  is  used  to  fit  the  points  in  the  data 
set.  This  minimizes  the  sum  of  the  squares  of  the  differences  between 
the  polynomial  and  the  data-set  values.  Two  solution  methods  are  used, 
"normal"  least  squares  and  orthogonal  polynomials.  Using  "normal"  least 
squares  assumes  a  solution  of  the  form 

Y  =  A0  +  ApX  +  A2X2  +  . . .  +  AnXn  .  (A-l) 

The  least-squares  procedure  leads  to  a  system  of  n  4-  1  equations  which 
are  solved  for  the  values  of  A  by  the  use  of  a  double  precision  matrix 
inversion  routine.  This  method  begins  to  lose  its  stability  above 
degree  11,  and  an  alternate  solution  technique  must  be  used. 

Orthogonal  polynomials  assumes  a  solution  of  the  form 

Y  =  ao  +  alQl  +  a2^2  +  • • •  +  anQn  5  (A-2) 
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TABLE  A- I 


Test  No.* 
258-10 
258-11 
258-12 
258-13 
258-14 
258-15 
258-16 
258-17 
258-18 
258-20 
258-21 
258-23 
258-24 
258-25 
258-26 
258-27 


First  and  Second  Order  Equations  of  Drag  Coefficient 
(Cd)  Versus  Velocity 


First-Order  Equation 
CD  =  3.469  -  0.0295V 
CD  =  1.837  -  0.00978V 
CD  =  3.115  -  0.0248V 
CD  =  2.382  -  0.0156V 
CD  =  1.917  -  0.0104V 
CQ  =  3.267  -  0.025V 
CD  =  1.613  -  0.0070V 
CD  =  1.503  -  0.00585V 
CD  =  2.421  -  0.0132V 
CD  =  2.982  -  0.0227V 
CD  =  4.6159  -  0.04104V 
CD  =  3.020  -  0.0214V 
CD  =  2.352  -  0.01586V 
CD  =  2.674  -  0.01844V 
CD  =  2.811  -  0.01968V 
CD  =  3.584  -  0.0292V 


_ Second-Order  Equation _ 

CD  =  12.58  -  0.257V  +  0.0014V2 
CD  =  6.76  -  0.124V  +  0.00067V2 
CD  =7.83  -  0.134V  +  0.00063V2 
CD  =  -2.95  -  0.106V  +  0.000695V2 
CD  =  3.03  -  0.035V  +  0.000135V2 
CD  =  19.17  -  0.39V  +  0.00209V2 
CD  =  -0.0263  +  0.039V  -  0.00031V2 
CD  =  -0.0599  +  0.039V  -  0.00031V2 
CD  =  -3.23  +  0.016V  -  0.00136V2 
CD  =  3.673  -  0.042V  +  0.000135V2 
CD  =  1.799  +  0.0442V  -  0.000629V2 
CD  =  2.309  -  0.00037V  -  0.000150V2 
CD  =  1.404  +  0.00952V  -  0.000164V2 
CD  =  7.089  -  0.1179V  +  0.0005565V2 
CD  =  -6.240  +  0.2031V  -  0.00137V2 
CD  =  -0.1097  -  0.2184V  -  0.00121V2 


“No  velocity  data  obtained 


for  Tests  258-19 


and  258-22. 
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Figures  A-l  through  A- 3 2 

Curve  fits  of  drag  coefficient 
versus  velocity  and 
error  versus  velocity 


NOTE:  Velocity  is  in  feet  per  second. 
Error  is  in  percent . 
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where  the  Q's  are  each  polynomials  which  are  orthogonal  to  each  other. 
They  are  defined  by  a  three-term  recurrence  relation: 

=  XQk(X)  -  «k+1Qk(X)  -  ^k+1Qk_! (X)  .  (A-3) 


E  XiQK(Xi) 
i _ 

E  Q|(Xi) 

1  (sums  over  data  points)  (A-4) 

SxiWWxi> 

1 

Thus,  the  solution  in  this  case  yields  not  only  a' s  but  also  the  oe's 
and  /?fs.  This  method  is  stable  up  through  degree  50  (maximum  allowed 
by  the  program) . 

Any  zone,  except  the  first,  whose  fit  is  made  using  the  "normal" 
least-squares  method  may  be  constrained  by  forcing  the  value  at  a  given 
point,  and  the  derivative  at  that  point,  to  be  the  values  as  determined 
by  the  curve  in  the  previous  zone.  The  mechanics  for  doing  this  involve 
the  use  of  the  Lagrange  multipliers.  Two  functions,  g  and  h,  are  de¬ 
fined  as 

g  =  A0  +  AxXq  +  A£Xq  +  . . .  +  AnXn  -  Y0  =  0  (A-5) 

h  =  AL  +  2A2Xq  +  3A3Xq  +  ...  +  tiA^Xq-1  -  Yq  =  0  ,  (A-6) 

where  Yq  is  the  value  at  Xq  defined  by  polynomial  over  the  previous 
zone,  Yq  is  the  derivative  defined  in  same  manner. 

The  function  to  be  minimized  is  then  assumed  to  be 

G  =  f  +  2Ag  +  2fih  ,  (A-7) 

where  A  and  n  are  the  Lagrange  multipliers.  The  system  of  equations  in 
Ak,  A  >  a^d  U  are  generated  in  the  normal  manner.  These,  together  with 


With  Q0 (X)  =  1, 


Qi(x) 


a 


K+l 


£ 


K+l 
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Eqs .  (A-6)  and  (A-7)  are  solved,  and  the  resulting  A's  are  the  coeffi¬ 
cients  to  the  constrained  polynomial.  Although  A  and  n  are  solved  for, 
they  are  used  only  in  the  formulation  and  not  in  the  evaluation  of  the 
resultant  polynomial.  It  is  noted  that  the  previous  zone  may  have  been 
fitted  using  either  "normal"  or  orthogonal  techniques . 

Zone  determination  may  be  assisted  by  Figure  A-33 .  Although  it 
is  only  a  rough  indication,  acceptable  results  have  been  obtained  using 
this  scheme  on  all  data  sets  tested. 


The  error  determination  is  defined  as 


E 


Yi  -  Y<Xi> 
Y(Xi) 


(A-8) 


where  Y(Xi)  is  the  data  value  and  Yj_  is  the  computed  value  using  the 
fitted  polynomial . 
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